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I. INTRODUCTION 

The transition from a confined phase with broken 
chiral symmetry to a deconfined chirally symmetric 
phase is an important subject for studies of finite 
temperature quantum chromodynamics (QCD). This 
transition is relevant for the evolution of the early uni- 
verse and reproduced in current heavy ion collision ex- 
periments. It can be investigated non-perturbatively 
using lattice QCD as long as the chemical potential 
for fermion number is small, ji/T < 1. A lot of ef- 
fort has been invested in lattice studies at zero chem- 
ical potential, for recent reviews see [H-Ql. Impres- 
sive progress has been reported very recently by sev- 
eral collaborations working with different fermion dis- 
cretization schemes (sl-fiot. In particular, lattice QCD 
with staggered fermions and physical quark masses 
does not predict a true phase transition but an an- 
alytic crossover in the limit of zero chemical poten- 
tial pl|. Similarly, results on the transition tem- 
perature and the equation of state have predomi- 
nantly been obtained from simulations with staggered 
fermions However, this fermion discretization 

is subject to an on-going debate and there is no formal 
proof that its continuum limit will reproduce the uni- 
versality class of QCD [l3] • It is therefore desirable to 
obtain independent results with other discretizations, 
in order to have some mutual control over systematic 
errors. 

Unfortunately, Wilson-type fermions (and even 
more so chiral fermion formulations) require higher 
computational costs. It is thus expedient to study 
the nature of the phase transition for various larger 
than physical quark masses and to extrapolate to the 
physical situation. Moreover, knowing global proper- 



ties of the phase transition as a function of the light 
quark masses constrains the enlarged phase diagram 
including the strange quark and non- vanishing chem- 
ical potential [ll]. An as yet unsettled crucial ques- 
tion in this context is the nature of the phase tran- 
sition in the two-flavor chiral limit. Most studies fa- 
vor a second order transition in the 0(4) universality 
class [T9l-[23but there are also claims for a first order 
transition~ j25l - [29j . Since in the continuum and chiral 
limits the transition is associated with the breaking of 
a global chiral symmetry, it is necessarily a true and 
non-analytic phase transition and one of these scenar- 
ios has to be realized [S^l, while an analytic crossover 
is ruled out. On the other hand, for moderate and in- 
termediate quark masses, the transition is an analytic 
crossover, before it turns into a first-order deconfine- 
ment transition for very heavy quarks. 

In this article we study the thermal transition with 
two degenerate flavors of maximally twisted mass 
fermions, which provide an C'(a)-improved Wilson 
fermion discretization, for a review see [3l|. As 
a first step we focus on the determination of the 
phase boundary, i. e. the pseudo-critical temperatures 
7c(m7r) using the Polyakov loop, the chiral condensate 
and the plaquette as observables. We do this for a set 
of pion masses, TOtt ~ 300 — 500 MeV, and attempt 
various extrapolations to the Nf = 2 chiral limit. Sim- 
ilar efforts were recently under way employing clover 
improved fermions [33 - [33 | . 

The following section serves to specify our simula- 
tion setup. In Sect ion Hill we introduce the observables 
and collect the pseudo-critical couplings from our sim- 
ulations. These results allow for an estimate of the size 
of the discretization errors present in our simulations. 
In Section [TV] we use these pseudo-critical points for 
an extrapolation to the chiral limit. We discuss pos- 
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sibilities and limitations in discerning the order of the 
chiral phase transition. Finally Section IVl gives some 
conclusions and an outlook. 



II. SIMULATION SETUP 

We consider QCD with a mass-degenerate dou- 
blet of twisted mass fermions, cf. the review by 
Shindler 31]. The gauge action is tree-level Symanzik 
improved while the fermion action is 

Sf[U,'iP,1P] ^J2x{x) (1 - kDw[U] + 2iKa/io75r^) xi^) 

X 

(1) 

The fermion fields are written in the twisted basis 
{X: x} which is commonly used for numerical simu- 
lations. It is connected to the basis of physical fields 
{■0,^} for the relevant case of maximal twist via 

-0 = -^(1 + i75'r^)x and V = + i75T^) ■ 

(2) 

The quark mass is determined by the hopping param- 
eter K, which parameterizes the untwisted bare quark 
mass component, 



K = (2amo -¥ 8r) , 



(3) 



and the twisted mass parameter /ip- The Wilson co- 
variant derivative is given by 



Y.{{r--i^)U^{x)ij{x + fL) 

(4) 



+ {r + ^^,)Ul(x-(i)i;(x--fi)) . 

In the weak coupling limit, /3 = G/^q — > oo, zero quark 
mass corresponds to k = 1/8, setting r — 1. For finite 
coupling this value of n gets corrections through mass 
renormalization. The overall renormalized quark mass 
M is composed of the twisted and untwisted masses 
as 



Z'L {mo - m„.) + ZIpI . 



(5) 



At maximal twist, the above fermion formulation is 
automatically 0(a)-improved, i. e. cutoff effects linear 
in the lattice spacing a are absent for non-zero physical 
observables. Maximal twist is achieved by tuning the 
hopping parameter to its critical value Kc, correspond- 
ing to mcr, where the untwisted theory would feature 
massless pious. The required knowledge of Kc{I3), as 
well as other input needed from zero temperature sim- 
ulations in order to set the scale, can be interpolated 
from data by the European Twisted Mass Collabo- 
ration (ETMC) [111. In Figs. [U and [2] we show our 
interpolations for /tc(/3) and the lattice spacing a(/3). 
Our numerical evaluation proceeds by an HMC algo- 
rithm [36^ within the publicly available code for QCD 
with twisted mass fermions [37]. 
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FIG. 1: Interpolation of the critical hopping parameter 
from ETMC data. 
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FIG. 2: Interpolation of the lattice spacing from ETMC 
data. 



Wilson fermions are well-known to feature unphys- 
ical phases for light quarks and coarse lattice spac- 
ings. Like the physical parameter space, these get ex- 
tended to a third direction because of the additional 
twisted mass parameter in the current formulation. 
In order to stay away from unphysical regions, knowl- 
edge of the bare parameter phase diagram is required, 
which we have mapped out earlier in a preparatory 
study [s^ . Status reports of our ongoing project have 
been given at the annual lattice conferences ]39[ [ioj . 

The temperature scale is set by the temporal lattice 
extent and the lattice spacing, T = \/{aNr). In order 
to locate the phase boundary between the hadronic re- 
gion and the quark giuon plasma, we perform scans in 
the lattice gauge coupling /3, which thus corresponds 
to a change in temperature of the lattice system. Ta- 
ble |T] gives the list of runs for different pion masses and 
the naming scheme that we have adopted for the sake 
of simplicity. To adjust the masses, ETMC provides 
parameters for NNLOxpt-formulae at their values of 
/3 e {3.8,3.9,4.05,4.2} which can be used to iden- 
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Run 


iV^ X 


Range 


(MeV) 




A12 


32^ X 12 


3.84 <I3< 3.99 


316(16) 


0.673(42) 


B12 


32^ X 12 


3.86 <I3< 4.35 


398(20) 


0.847(53) 


C12 


32^ X 12 


3.90 <I3< 4.07 


469(24) 


0.998(62) 


BIO 


32^ X 10 


3.76 <I3< 4.35 


398(20) 


0.847(53) 



TABLE I: List of scans in /?. See also table|V] 



tify the relation m7r(/io) at those couphngs. For our 
Nr = 12 scans we have relied on the one- loop scaling 
relation 



a/xo(/3) = C exp 



12/3o 



(6) 



with Pq = (11 - 2iV//3)/(47r)2 and fixing the free 
parameter C at one of the available couplings. We 
have found this relation to work sufficiently well to 
create lines of constant pion mass within the errors 
quoted in TablelH The run at Nr = 10 has a constant 
a/io = 0.006 in the /3-interval from 3.865 to 3.930 for 
which we likewise have the same pion mass within er- 
rors in our simulation range. For the other /3- values 
we have adapted the twisted mass according to a two- 
loop scaling relation similar to the one-loop formula 
shown above. The free parameter C has been adapted 
to produce ano = 0.006 at /3 = 3.88. 

A final comment concerns the explicit flavor sym- 
metry breaking due to the twisted mass term at finite 
values of the lattice spacing. This breaking has been 
investigated by the ETM collaboration for T = the- 
oretically [41)] and in simulations [HI . The outcome is 
that effects from flavor breaking - formally of 0{a^) 
- appear to be negligible in all quantities investigated 
so far but the neutral pion mass. For this reason we 
use the charged pion mass throughout the paper. As 
will be explained in Section ITVl for our scaling analysis 
we need to be close enough to the continuum in order 
to reproduce chiral symmetry, where flavor breaking 
should not play any role any longer. Comparison of 
two lattice spacings appears to justify this assump- 
tion. However, a third value of the lattice spacing is 
required in order to make these statements about the 
size of lattice artifacts more definite. 



III. THERMAL TRANSITION 
TEMPERATURE 

In order to locate the transition, we have used both 
pure gauge and fermionic observables. The gauge ob- 
servables are the plaquette 



P 



1 



I^eTr^^C/^,(x) 



(7) 



with 

U^,,ix) = U^{x)U,{x + fi)Ul{x + v)Ul{x) , (8) 
and the real part of the Polyakov loop 
1 1 



(9) 



X 0:4— 



The latter is of particular interest since it is the order 
parameter of the pure gauge deconfinement transition. 
Along with these observables, we look at their suscep- 
tibilities, 



Xo-Nl({0^)-{Of) 



(10) 



The renormalized (real part of the) Polyakov loop 
can be determined as [i^l 

(Re(L))^ = (Re(L)) exp {V(r^)/2T) , (11) 

where V{rQ) denotes the static quark- antiquark po- 
tential at the distance of the Sommer scale r = tq [43| 
to be determined at zero temperature. 

The chiral condensate (?/"/') represents the real or- 
der parameter of chiral symmetry breaking in the 
massless limit. An appropriate quantity to locate the 
chiral phase transition is the chiral susceptibility 



dm„ 



(12) 



Here, we consider only a part of that expression, the 
variance per configuration, 



VlT({{^^f)~{H,f) . (13) 



This quantity shows a peak associated with the chiral 
transition. Moreover, it is expected to dominate the 
signal of Xct, see e.g. [43 |. 
The pion norm 

Kl' = E (V^(^)^75T+V(a:) V^(0)i75T-V(0)\ (14) 

is interesting for twisted mass simulations because its 
definition is independent of the fermion basis. It is 
connected with the chiral condensate via 



(15) 



which has been proven for lattice twisted mass 
fermions in [4^ . We have used this relation as a check 
for (V'V')- 

At maximal twist the chiral condensate can be 
renormalized as follows (see the appendix in [46i | and 
references cited therein) 



(^V)^ = Zp((^V')+c(ffo)^). 



(16) 
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This immediately suggests the form of a subtracted 
condensate, which is completely standard. However, 
the subtracted condensate is no longer an order pa- 
rameter for the chiral transition. It is very easy to fix 
this problem by adding the zero temperature chiral 
condensate in the chiral limit. Thus, we introduce a 
(re)normalized condensate in terms of the ratio 

„_ ^ (?V')(r,Mo) - W)(o,Mo) + W)(o,o) 
W)(o,o) 

(17) 

where {^il?j{T, /io) means ('0V') to be evaluated at non- 
zero temperature and finite /zq. {'ip'ip){Q, /io) can be ob- 
tained from spline interpolations of T = (V'V') data 
in both the mass /xq and /3. Additionally, to determine 
{tpip){0, 0) one has to perform a chiral extrapolation of 
the T = (tptp) data at every /3- value one is interested 
in. We have used a linear extrapolation through three 
points at every /3. The data turned out to be com- 
patible with a linear fiQ dependence over the whole 
temperature range we consider here. For the T = 
data we were relying on results provided by the ETM 
collaboration. 

The fermionic observables have been determined us- 
ing the technique of noisy estimators, as in For 
|7rp we have calculated ten propagators per gauge con- 
figuration on Z{2) noise vectors. (V'V') was evaluated 
using 24 Gaussian volume source vectors for BIO, B12 
and C12, and 24 Z(2) volume source vectors for A12 
respectively. All propagators have been calculated on 
commodity graphics hardware using NVIDIA's CUDA 
programming language. The statistics accumulated 
for the various runs as well as the averages for the 
Polyakov loop and the chiral condensate are given at 
the end of the paper in TablelVl 

Quite generally, we find the signals for the transition 
to be quite smooth and noisy which presumably is re- 
lated to the fact that we are merely probing a very soft 
crossover in our range of pion masses. For a crossover 
there is no unique definition of a critical temperature 
as the physics changes smoothly and analytically be- 
tween the different regions. A pseudo-critical temper- 
ature associated with the transition behavior of indi- 
vidual observables is in general observable-dependent. 

In Figs.[3]-E]we show our data for a— in accordance 

with Eq. (IT51) and the susceptibility of the real part 
of the Polyakov loop (Eq. ©). We quite clearly see 
maxima for in all cases, whereas for the Polyakov 
loop susceptibility we find only an onset of certain 
shoulders for the ensembles A12, B12 and BIO. At 
the higher pion mass case C12, where we restricted 
ourselves to smaller statistics, there seems to appear 
a maximum also for the Polyakov susceptibility. 

In order to estimate the pseudo-critical /3c for chiral 

transition we have modeled the data for a?- , with a 

ip'ip 



Run iV^ 


13c 


Tx (MeV) 




A12 12 


3.89(3) 


202(7) 


0.437(18) 


B12 12 


3.93(2) 


217(5) 


0.473(10) 


C12 12 


3.97(3) 


229(5) 


0.500(14) 


TABLE II: List of pseudo-critical points for the chiral tran- 


sition T^. 








Run iV^ 




Tdeconf (MeV) 


T'oTdcconf 


B12 12 


4.027(14) 


249(5) 


0.546(13) 


C12 12 


4.050(15) 


258(5) 


0.565(14) 



TABLE III: List of pseudo-critical points for the decon- 
finement transition Tdeconf. 



Gaussian 

c + aexp(^-^^-^^ . (18) 

The results for the corresponding pseudo-critical 
chiral transition temperature are collected in Ta- 

bienn 

In Fig. [7] we show the renormalized chiral conden- 
sate ratio R(^^-^ and the renormalized Polyakov loop 
(Re(L)) for the ensembles B12 and BIO. The large er- 
ror bars for the T-values in case of the BIO ensemble 
reflect the uncertainty in the scale setting. 

By determining the inflection point of the renormal- 
ized Polyakov loop (Re(L)) we were able to estimate 
the deconfinement temperatures Tdoconf for the ensem- 
bles B12 and C12, see Table|llll 

We cleary see that Tdeconf > ?x f*-"^' both higher 
pion masses. This corresponds to the observation re- 
ported in [42]. 

From weak coupling analyses (valid at high temper- 
ature) it is known that the leading order a^-scaling 
towards the continuum limit might not set in before 
^ 16 [48]. Therefore, discretization effects as a 
major source of systematical errors need to be thor- 
oughly checked. Since the runs BIO and B12 share a 
common pion mass and differ only by Nr , they can be 
used in order to assess the magnitude of cutoff effects. 
As can be seen from Fig. [6] the quality of <j'^^ for the 

BIO ensemble is not yet precise enough to allow for a 
Gaussian fit. The available data however suggests a 
maximum at around /3 ~ 3.82 which corresponds to 
a temperature T ^ 218 MeV and agrees with T-^ at 
Nr = 12. Moreover the renormalized Polyakov loop 
and the renormalized chiral condensate (Fig. [7]) agree 
within errors for BIO and B12 indicating small cutoff 
effects. 
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FIG. 3: a-|^ (left) and susceptibility of Rc(L) (right), both for run A12. 
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FIG. 4: a^^ (left) and susceptibility of Re(L) (right), both for run B12. 
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FIG. 5: a-^^ (left) and susceptibility of Re(L) (right), both for run C12. 
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FIG. 6: cr|^ (left) and susceptibility of Re(L) (right), both for run BIO. 
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FIG. 7: Ratio R(:^^) acc. to Eq. (fT7)l (left) and renormalized Polyakov loop (Re(L))^ (right), both for runs B12 and 
BIO. 



IV. TOWARDS THE CHIRAL LIMIT 



As indicated in the Introduction, the main inter- 
est in the Nf = 2 thermal transition lies in its chiral 
limit, for which one would like to unequivocally de- 
termine the order of the phase transition. The chiral 
condensate (ipip) then is an order parameter corre- 
sponding to the magnetization in an appropriate spin 
model of the same universality class. Finite quark 
(and therefore pion) masses break the chiral symme- 
try explicitly, thus corresponding to an external field. 
Provided the Nf — 2 chiral limit features a second 
order transition and belongs to the 0(4) universal- 
ity class, one may extrapolate finite mass simulations 
using universal scaling relations, which hold within 
some scaling region around the critical phase transi- 
tion [l^ [2l| . A priori it is not known how far into the 
massive region scaling extends, i.e. one can merely 
test consistency of the data with scaling. A further 
difficulty is that chiral symmetry is broken explicitly 



for Wilson fermions at finite lattice spacing, even in 
the massless case. Any universal behavior for these 
types of fermions thus corresponds to continuum scal- 
ing, which can only be observed once discretization 
errors are sufficiently small. Finally, the scaling rela- 
tions we employ here are valid in the thermodynamic 
limit. Dedicated finite-size scaling analyses are re- 
quired to establish the appropriate lattice sizes, but 
this is beyond the scope of the present study. Again, 
we assume our lattices to be sufficiently large and test 
for consistency with scaling. 

We begin by attempting a fit of T^{mT^) to the scal- 
ing form [2l|,|4| 

r^(™.)=T^(0) + ^-m^/(^^) , (19) 

where we have dressed the critical exponent jS with 
a tilde in order to distinguish it from the lattice cou- 
pling. The "external field" in this case is the quark 
mass specified by the mass parameter a^Q, which in 
turn is connected to the pion mass in LOxpt via 
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^ /io- Thus, it is important to keep the pion 
mass smaU for two reasons, the vahdity of both the 
scahng window and the LO of xpt- While there is 
good reason to expect that our pion masses are suffi- 
ciently small for the latter [s^ , the size of the scaling 
region remains unknown at present. Unfortunately, 
we do not have sufficiently many data points or suf- 
ficiently small errors in order to determine the expo- 
nents, but fix the exponents and fit A and Tc{0) only. 
For 0(4) we have 2/{f35) = f .08 and the resulting ex- 
trapolation is shown in Fig. [51 giving a chiral critical 
temperature T^{mTr = 0) = 152(26) MeV. It is now 
interesting to ask whether 0(4) scaling can be dis- 
criminated from other behavior. As discussed earlier, 
the alternative scenario is a first order phase transi- 
tion in the chiral limit. Often in the literature the 
same scaling relation is tested by merely c hang ing to 
"first order exponents" (2/(/35) = 2) 0, Hof. Do- 
ing so leads to an extrapolation with somewhat larger 
T^{m^ = 0) = 182(14) MeV. However, it is unclear 
to us whether the scaling relation is applicable in this 
case. Firstly, for a first order phase transition there 
is no diverging correlation length. Approaching in 
the infinite-volume limit from above and below pro- 
ceeds in different phases, with finite correlation length 
in each. Hence, there is no scaling and no universality 
in the sense of second order transitions (in particular 
(3 = and S = oo separately). The "critical expo- 
nents" usually associated with first order transitions 
specify the approach of the thermodynamic limit in 
finite-size scaling analyses, but do not apply to the 
relation (jl9p in the thermodynamic limit (for a de- 
tailed discussion of scaling for first order phase tran- 
sitions, see [HI). Secondly, if the chiral limit indeed 
features a first order phase transition, it will weaken 
with finite quark masses until it vanishes in a Z{2) 
critical endpoint. Fig. [HI shows the two possible sce- 
narios. However, this means that coming from the 
crossover region at larger quark masses, an extrapola- 
tion to the chiral limit is never exact, as it would pass 
through a singularity at the critical point. Rather, the 
approach of this singularity will again be characterized 
by scaling, this time in the Z{2) universality class. In 
this case we may use again the relation ([19]), but with 
a finite critical pion mass marking the critical point, 
toJ — ^ (m^ — J,). We have attempted such extrap- 
olations also. Our data are not sufficient to constrain 
TOjr.c- Therefore, Fig. [8] shows two extrapolations, one 
with m7r,c ~ and another with m^r^c ~ 200 MeV. As 
the figure illustrates, our extrapolations alone cannot 
yet discriminate between the first order and second or- 
der scenarios. This would require drastically smaller 
pion masses, lower than the physical value even. Nev- 
ertheless, utilising knowledge about Tc from other sim- 
ulations we still obtain a tendency. The fit assuming 
a first order scenario leads to a critical temperature 
which is somewhat larger than expected from other 
investigations ll]. Of course, those extrapolations are 
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FIG. 9: Illustration of possible scenarios for the Nf = 2 
chiral limit. 



likewise valid only in the 0(4) scenario, so again this 
is merely a consistency test. 

For a fixed Nt, assumed to be large enough so as 
to be sufficiently close to the continuum, it is also 
possible to obtain the chiral critical /3 by means of the 
scaling relation [2l|, |49j 



I3c{h) = /^chiral + B ■ h'/<^^'^ 



2a/io (20) 



with 1/{(3S) = 0.537 corresponding to 0(4) exponents. 
For = 12 our estimates for (3c are shown in Fig.[TUl 
and can be extrapolated in this manner. 

Consistent fits have been found taking all three 
points from A12 to C12 into account. The result for 
the critical chiral /3-value is 



/3chiral(iVr = 12) = 3.73(9) . 



(21) 



We have carried out the same fit but with the two 
lower pion mass values (A12 and B12) only. It ended 
up with the same value. This result corresponds to 
T^{m^ = 0) « 152(26) MeV where the error results 
from he scale setting. This number is in accord with 
our fits for T-^{mT^) for a second order transition in the 
chiral limit. Note however, that the lattice spacing 
necessary to set the scale stems from an extrapolation 
to smaller values of /3 than available from ETMC. This 
is refiected in the large uncertainty assigned to the 
temperatures. 

Next, the scaling of the magnetic equation of state 
can be investigated, Fig. llli where we follow previous 
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FIG. 11: Scaling for the bare {^ip) data at A'^^ = 12 as 
function of the scaling variable with modelling of scaling 
violations. The fit shown is for the combined A12 and B12 
data (fit No 10 in Table HVl . 
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with 



of /3chirai wlth increasing mass. Indeed, B12 yields 
/3chirai = 3.40(5), which would correspond to a much 
too low critical temperature below 100 MeV, while 
C12 gives even smaller values with larger . Thus, 
the fit seems to account for scaling violations due to 
large mass by decreasing /^chirai- 

However, scaling violations due to the quark mass 
can be taken into account by an ansatz including cor- 
rections (2^ . 

(V^V) = h^'^cf{dT/h^l'^~^^^) + atTh + bih + b:ih^ + ... . 

(24) 

We have fitted our data in numerous ways by tak- 
ing into account one, two or even three violation 
terms. Joint fits to the A12 -|- B12 ensembles are 
feasible in all three combinations, giving a /3chirai the 
more consistent with the previous determination the 
more violation terms are included (see TablelTV)) . In 
Fig-El we show a combined fit to A12 and B12 fix- 
ing /3chirai = 3.73 from our independent determination 
with x^/dof = 0.63. Note that these fits with the two 
lower order violation terms are not able to include the 
C12 data with the requirement of a reasonable value 
of x^/dof. However, if we include the next higher vio- 
lation term h^h^ in the combined fit to A12, B12 and 
C12 we obtain an acceptable x^/ dof — 1.8, see the 
the last line of Table lTVl We observe that in this case 
the fit even prefers a value for /^chirai compatible with 
the one from the analysis based on eq. ([20l) . 

Since we are in a range of the scaling variable 
r/Zi^/'^^'^) where the scaling function is rather flat, 
judgement on whether there are additional violations 
of the 0(4) behavior or not is difficult. Repeating 
this exercise for the first order scenario with end- 
point does not give further insight as the combinations 
of exponents are very close, l/(/3(5) = 0.537,0.638 
and 1/5 = 0.21,0.20 for 0(4) and Z(2), respectively. 
Therefore, our data are consistent with the 0(4) sce- 
nario, but do not rule out the possibility of the first 
order case. This would require drastically smaller pion 
masses combined with finite-size studies, as the win- 
dow for chiral scaling appears to set in for m.,^ <^ 300 
MeV. 



r = /3 - /3chiral . (23) 

The functional form of the scaling function / for 
the 0(4) case is known [5^ [5^. Since we do not 
know the correct normalization for r and h with re- 
spect to QCD, we are left with two free parameters, 
c and d, that have to be fitted. We perform the 
fits in the /3-intervals from /3 = 3.83(3.85,3.89) to 
P = 3.97(4.03,4.04) for A12 (B12, C12), respectively. 
The fit results are collected in Table lTVl 

A fit for the line A12 works quite well, x^/dof = 
0.43, but giving /3chirai = 3.57(4), which is smaller 
than the value estimated above by applying eq. (|20|) . 
In general we observe an increase of x^ and a decrease 



V. CONCLUSIONS 

We have presented a (revised) first investigation 
of the two-flavor thermal QCD transition with maxi- 
mally twisted mass fermions. Our results are compat- 
ible with existing work although, of course , stagg ered 
investigations are much more advanced [T2l . ll3l[l5[ , [53 |. 
The quality of our signals is comparable to recent re- 
sults with clover improved Wilson fermions [H, [s^l • 

For three pion masses in the range 300 MeV < 
m^r < 500 MeV we have determined pseudo-critical 
temperatures for the crossover from the hadronic 
regime to the quark gluon plasma. The pseudo-critical 
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No 


Data 


/-'cniral 




d 


at 








1 


A12 


3.57(4) 


0.14(2) 


0.367(7) 











0.43 


2 


B12 


3.40(5) 


0.22(4) 


0.36(2) 











0.64 


3 


C12 


3.12(2) 


0.42(3) 


0.39(2) 











2.42 


4 


A12 + B12 


3.368(6) 


0.257(6) 


0.383(5) 











3.31 


5 


A12 + B12 


3.48(2) 


0.225(6) 


0.48(2) 


0.7(1) 








2.2 


6 


A12 + B12 


3.57(2) 


0.152(7) 


0.53(2) 





0.90(6) 





1.75 


7 


A12 + B12 


3.82(4) 


0.028(9) 


1.1(2) 


-2.2(2) 


2.49(8) 





0.42 


8 


A12 + B12 


3.73 


0.1279(8) 


0.825(8) 


4.01(4) 








76 


9 


A12 + B12 


3.73 


0.0759(7) 


0.81(2) 





1.61(2) 





7.2 


10 


A12 + B12 


3.73 


0.053(2) 


0.74(2) 


-1.8(2) 


2.23(6) 





0.63 


11 


A12 + B12 + C12 


3.76(2) 


0.047(6) 


0.83(6) 


-1.5(2) 


2.20(6) 


50(11) 


1.8 



TABLE IV: Fit results based on eq. (|24[) for several combinations of our data sets and fit parameters. Numbers in bold 
face have been fixed before fitting. The fit shown in Fig. l 1 1 1 corresponds to line No 10. 



temperatures - extracted for the two higher mass val- 
ues - from observables related to chiral and deconfine- 
ment transitions, respectively, turned out to be differ- 
ent. Discretization effects in appeared to be small 
for our lattice spacings, a < 0.09 fm. 

We have restricted ourselves to pion masses < 
500 MeV in order to assure the validity of LOxpt as 
well as the scaling forms in order to extrapolate to the 
chiral limit. Assuming the scaling forms appropriate 
for different universality classes, such extrapolations 
gave critical temperatures in the range Tc ~ 140 — 200 
MeV consistent with other studies. However, detailed 
fitting analyses demonstrated that the second order 
0(4) scaling regime is not yet reached. Scaling viola- 
tions could be accomodated by leading order correc- 
tions due to finite-mass effects up to ttItt ^ 400 MeV, 
while heavier masses violate even those corrections. 
By including higher order violation effects reasonable 
fits could be achieved with /3c values consistent with 
the other determinations. 

We find that truly distinguishing between the dif- 



ferent universality classes and thus ruling out a first 
order scenario will require much smaller pion masses, 
m^r < m^^^^ as well as finite-size scaling analyses. We 
hope to address these issues in future investigations. 
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TABLE V: Statistics for gauge observables from our simulations as well as expectation values of Re(I/) and {tptp). Note 
that the trajectory length differs between the runs. On the apeNEXT (BIO, C12 except for /3 = 4.06) t = 0.5, on the 
HLRN (A12, B12, /3 = 4.06 of C12) r = 1. 
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